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Abstract 

Simple derivation of Spitzer's asymptotic law for Brownian windings 
(Trans Am. Math. Soc. 87 ,187 (1958)) is presented along with its generalizations. These 
include the cases of planar Brownian walks interacting with a single puncture and Brownian 
walks on a single truncated cone with variable conical angle interacting with the truncated 
conical tip. Such situations are typical in the theories of quantum Hall effect and 2+1 quan- 
tum gravity respectively. They also have some applications in polymer physics. Extension 
of these results to the the multiple punctured case is also briefly discussed. It is technically 
associated with some results known in the context of string and conformal field theories and 
theories of quantum chaos. 
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In 1958 ,based on earlier work by Paul Levy [1] ,Spitzer [2] had obtained the asymtotic 
probability distribution P(x) for the winding angle 9 for the planar Brownian motion. If 
z(t)=x(t)+iy(t) is 2 -dimensional Wiener process,then it is of interest to study distributions 
of \z(t)\ and arg z(t)=9(t) .For large times t , Spitzer obtained his famous Cauchy-type 
distribution for 9{t) given by 

29 1 1 

P(x = - — )dx = -dx . (1) 

y lnr irl + x 2 w 

Eq.(l) is obtained under the assumption that the random walker begins his travel at some 
point z in z-plane other than the origin , i.e. x=y=0 . Then,the angle 9 in measured 
with respect to the line which joins z with the origin. This problem is of interest in poly- 
mer physics [3-5] since it represents the benchmark problem for study of entanglements .In 
mathematical literature the same phenomenon is described in terms of recurrence and tran- 
sience. For example, it is well known [6] ,that 1- and 2- dimensional Brownian motions are 
recurrent (that is the random walk visits time and again its starting point) while 3- dimen- 
sional motion is transient (that is there is a nonzero probability that the walker will not 
return to the origin). For the random walk on once- punctured z-plane, Lyons and McKean 
[7] had demonstrated that the walk is recurrent while for the twice- punctured plane the 
walk is transient. In the language of polymers this means, that the polymer lying in z-plane 
will not be entangled with another polymer, placed perpendicular to this plane , while it 
will become entangled if there are at least 2 polymers which intersect z-plane at 2 distinct 
points. The planarity of the above problem is actually not to essential as it was explained 
in Ref.[5]. In the case of quantum mechanics the once puctured plane problem is directly 
associated with the Aharonov-Bohm (A-B) effect [8]. The A-B effect in the presence of 2 
punctures was studied in Ref.[9].The methods of Ref.[9] cannot be generalized to the case of 
more than two punctures and do not provide any information about the recurrence and/or 
transience. At the same time, the methods used by McKean and Lyons [7] can be used for 
the case of more than 2 punctures but are not widely known in physics literature. They 
had been recently mentioned in Ref. [5] in connection with some topological problems arising 
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in polymer physics. In physics literature random walks on multiply punctured plane were 
extensively studied in connection with problems related to the quantum Hall effect (QHE) 
and anyonic superconductivity [8, 10] while in mathematics literature the same problem was 
recently extensively studied by Pitman and Yor [11, 12]. To our knowlege,no attempt had 
been made to establish connections between these two formalisms (see also Ref.[7]) . In this 
brief note we would like to make the first step towards this comprehensive goal. 

Let us begin with the well-known expression for the distribution function for the planar 
random walk given by 

^i,r 2 ;t) = ^exp{- (ri ~ t r2)2 }, (2) 

where ri(r 2 ) = {xi(x2), yi (2/2)} -With respect to the origin in z-plane the polar system of 

coordinates can be used. In this system of coordinates Eq.(2) can be rewritten as 

1 2,2 00 

G(n,r 2 ,A6l;t)=— exp{— 4±-^} £ e imAe I m (z) , (3) 

,where A6 = 9 l — 62, z = 2rir 2 /t and I m (z) = I- m (z) is the modified Bessel's function. The 
above distribution function can be used either for study of the radial or the angular dis- 
tributions or both.Suppose,we are interested in the angular distribution function only (in 
view of Eq.(l) ). Then, using Eq.(3) ,it is convenient to introduce the normalized distribution 
function defined according to the following prescription: 
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The Fourier transform of such defined distribution function can now be obtained in a stan- 
dard way as 

f(z, a) = [ dA6 e— Ad f(z, A9) = ^£2 . (5 ) 

—00 

Let us now choose r 2 = fy/i + r\. This choice is motivated by known scaling properties 
of Brownian motion[13]. Then , for large t one obtains :z ~ 1r\rj\ft .For fixed f and r\ 
and t— > 00 one surely expects z — > 0. This observation allows us to use known asymptotic 



expansion for I\ a \{z) for small z's with the result for f(z,a) ( valid for small z's or large 
t's ): 

\ot\ 

f{z,a)^exp{ — 2"lnt}- (6) 

The inverse Fourier transform of Eq.(6) leads us to the result given by Eq.(l),i.e. f{z, A9) = 
P(x) where x=2A6>/lnt. 

Thus obtained result can be easily generalized now. For example, instead of considering 
random walks on the flat once punctured plane we can consider the same problem on the 
surface of a cone.This type of problem is of interest in connection with the study of 2+1 
quantum gravity [14,15]. It is also known [15,16] ,that the above conical problem is equivalent 
to the planar random walk problem in the wedge (the conical angle is simply related to that 
of the wedge). In the most general case our walk may be allowed to interact with the edges 
of the wedge [2] or, in the case of the wedge angle equal to 2n ,with the puncture located at 
the origin. The analogue of the distribution function ,Eq.(3), is known to be [16,17]: 

G(n,r 2 ,A9 ; t)=l expj-^} £ e^^t I, |m+ ,,(z) . (7) 

P l Zl m=-oo P 

For (3 = 2n and 5 = Eq.(7) reduces to Eq.(3) as required. The wedge angle (5 lies between 
and 2n while the statistics- changing parameter 5 is responsible for the polymer puncture 
interactions as it is explained in Refs.[5,18] or for the interaction with the flux tube if the 
magnetic language is being used [8]. 
By analogy with Eq.(4),we obtain, 
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A*) = j^j-: E e 2 ™-I 2flm+5{ (z) . (8) 
Upon Fourier transforming this expression we obtain, 

l2 f\s\( z ) 

Repeating the same chain of arguments which had lead us to Eq.(6) we obtain now, 



fp(z,a) ~ exp{-i( 



27T r 

a + J 5 



j mint}. (10) 



To perform the inverse Fourier transform of Eq.(lO) is nontrivial. Indeed, we have 
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f$(M,t) = ± J daexpiiaM -^(Ja 2 + (jS) -j\5\)\nt}. (11) 
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The integrals of this type are known in the context of quantum field theory [19] and had 
been also used recently in polymer physics problems [20]. By introducing new variable a = 
sinhy9 into Eq.(ll) it is transformed into 
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fp(A9,t) = — exp{— \5\ hit} / dtp cosh tp exp{ — - [cosh tp} — \5\ hit + iA6— \5\ sinh</?} . 
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(12) 

The exponent inside of the integral in Eq.(12) can be transformed as follows: 

\ 27r 2ix 

— -[coshy?]— \5\ lnt + iA9— \5\ smhtp = — \J a 2 + w 2 cosh((/? + tp ) , 
2 p p 

where a — | \5\ hit and uo = \5\ so that cosh(/9 = ~j^=p. an d sinh(^ = ^2^2 -Use of 

these results in Eq.(12) allows us to rewrite it in the equivalent form , 

. 5f\5\\nt vr J 

fg(A6,t) = — . exp{— hit} / dtpcoshtpexp[— V a 2 + u 2 cosh tp] 

flva 2 + u 2 p J 

= Q^K^VtfTu*) , (13) 

where Ki(x) is the modified Bessel's function with known asymptotic expansions: K\(x) — - 
for x and Ki(x) ~ \f^ e ~ x f° r x ~^ 00 • Using these expansions ,the following asymptotic 
results for the distribution function f^(A9,t) are obtained : 

1.5 — > and /3 is fixed and nonzero. In this case we recover Spitzer's law , Eq.(l) as 
required. 

2.5 — > 00 but AO is finite. Then,we obtain (x= 2 ^) , 

(3\nt 



x . . 1 / 5 dx r 
tnixjdx ~ -W — hit 5-exp|- 



x 2 } . (14) 
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Obtained result is in complete accord (up to numerical prefactors) with that obtained in 
Ref. [21] using methods of conformal field theory and, more recently, in Ref.[22] ,by direct 
numerical simulations. 
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3.5 and (3 are finite but AO — > oo and A# 3> lnt .Then, we obtain , 



f S p(x)dx ~ ^^lntexp{ 







lnt} exp{— |x| 
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InO^r- (15) 



2 



This result is in formal agreement with that obtained by Rudnick and Hu [4] and Desbois 
[23]. In Ref.[4] only the leading exponential factor was obtained while the result, Eq.( 15), 
differs from that obtained in Ref.[23] by an extra time- dependent factor .No connections 
with conical singularities or QHE problems were made in either Ref.[4] or Ref.[23]. 

With the results just obtained ,we are ready now to make connections with works of 
Ito and McKean [24] and Lyons and McKean[7] (see also Ref.[25]).We begin by stating the 
result proved by Levy[l] and refined by others [13]. Let z(t) be some planar Brownian 
motion which started ,say,at z =0 (at t = 0). Then , the motion obtained with help of 
some analytic function f(z) is also Brownian. This is equivalent to saying that the planar 
brownian motion is conformally invariant. Let us illustrate this fact in the example of once 
punctured plane R 2 — 0. 

The diffusion equation (in dimensionless units) on R 2 — can be written in a usual form 

as 

dt 4 \dx 2 dy 2 ) 1, 1 } 

where the factor | is chosen for convenience.lt is useful to rewrite Eq.(16) in terms of complex 
variables. Simple calculation produces 

v=*L. (it) 

dt dzdz V } 

Since this equation is not defined for z = z = 0,we introduce new complex variable w through 
z=exp w .Unlike z, our new variable w = u + iv is defined for the entire complex w— plane. In 
terms of w-variable Eq.(17) can be rewritten as 

Obtained equation describes Brownian motion on a simply- connected covering space M 
which is Riemann surface for the logarithmic function. The result just obtained coincides with 



that discussed in the book by Ito and McKean,e.g.see page 280 of ref.[24].Ito and McKean 
argue (without explicit demonstration) that , based on the results of P.Levy ,the diffusion 
Eq.(18) can be converted into standard form given by Eq.(16) by replacing Brownian time t 
by another (actually random) time T which is properly chosen. Let us demonstrate how this 
can actually be done. To this purpose, let us consider the Langevin-type equation written in 
the form of Ito : 

dx(s) = a(x(s))ds + o-(x(s))dw(s) . (19) 
The corresponding backward Kolmogorov-Fokker-Planck equation can be written now as 
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—P(x, s I so, so) = -a—P + -a 2 — 2 P . (20) 



ds ' ' dx 2 dx 

Introduce "new" time T(s) according to equation 

T(s) 



S = 





/ dtg(x(t)) . (21) 



Therefore, ds = g(x(T))<iT.Substitution of this result into Eq.(19) and use of Ito stochastic 
calculus [13] produces, 



dx(T) = a(x(T))g(x(T))dT + yj g(x(T))c(x(T))dw(T) . (22) 



New diffusion coefficient can be selected now ( in view of Eq.(20))as 

1 



D new = ^ga 2 . (23) 



Looking at Eq.(21) and selecting 



T(s) 

s = J dte 2u ® (24) 
o 



produces , in view of Eq.(23) , D new = | . Eq.(24) is in complete agreement with the result of 
Ito and McKean [21] where it was presented without derivation . Surely, T(s) is random time 
so that solution of Eq.(16) will depend upon random time. This is somewhat inconvenient. To 



correct this inconvenience it is useful to think about calculable observables.In view of Eq.(5) 
,we are interested in the Fourier-transformed angular distribution function f(z, a) .It can 
be shown [24], that this task is equivalent to finding averages of the type 



where (...) denotes the averaging with help of Gaussian-like propagator for the free "parti- 
cle" . Eq.(25) is a special case of famous Feynman-Kac formula[26] . The transition from first 
to second average in Eq.(25) is associated with the Cameron-Martin-Girsanov-type formula 
for changes of variables inside of path integrals[13,27].Eq.(25) is associated with Schrodinger- 
like differential equation. In view of Eq.s(21) and (24) ,it can be easily demonstrated , using 
standard methods of quantum mechanics [28] , that such obtained Shrodinger-like equation 
is just that for the modified Bessel functions , e.g. I m (x),etc .Its solution , indeed, leads us 
to Eq.(5). 

Based on these results ,immediate generalizations are possible. For instance, if we would 
be looking for the distribution function f(z, a) for 2- (or more) punctured plane,we would 
need the uniformizing function which will help us to write the diffusion equation on the 
simply connected universal covering surface M .This surface is known to be the Riemann 
surface for the punctured torus [25]. Classical motion on such surface is chaotic and quantum 
description of such motion was considered in Ref. [29] where no attempts were made to relate 
it to Spitzer's results. Already for two punctures study of transience and recurrence is highly 
nontrivial [7,25] . For two or more punctures the task of finding the uniformizing function can 
be reduced to that of finding all solutions of the corresponding Fuchsian-type equations [30] 
which are identical in the form to that known in string and conformal field theories for the 
correlation functions[31].It remains a challenging problem to extend Spitzer-like results to 
multiply connected surfaces,to connect these results with those of Pitman and Yor[ll,12] 
and with those known in the theories of QHE [8,10]. 

Note added in proof. When this work was completed, two additional recent results 
came to our attention. In Ref. [32] reader can find more up to date(as compared to Pitman 
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(25) 
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an Yor's) references on Brownian windings while Ref.[33] contains detailed applications of 
the above ideas to QHE 
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